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Abstract 


Initially spherical blast waves are systematically distorted due to persistent 
latitudinal solar wind structure. This is to be distinguished from the non-systematic 
(random) distortion due to varying structure in longitude which introduces an 30^ 
average deflection of shock normals from radial. The systematic latitudinal effect 
should be at least 25° at mid-latitudes, and observable in the 30° noise level with 
14 or more shocks for statistics. The observed occurrence rate of shocks during 
solar maximum is sufficient to detect the effect. Corotating shocks should become 
deter:tabie between 1 and 5 AU. identification could be a problem because of the 
30° noise level becoming greater beyond 1 AU. However, the three-dimensional 
geometry of corotating shocks show a strong latitudinal structure which can be used 
in an out-of-the-ecliptic mission for a statistical identification based on shock 
occurrence rates. 
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Introduction 


Most of the interplanetary shock waves observed with 1 AU of the sun 
originate from some short lived solar event, such as a solar flare, and then 
propagate out as a more-or*less spherical shock wove until they leove the solar 
system. Beyond 1 AU another class of interplanetary shock wave becomes common— 
the cortating shock pair formed by the interaction of long lived solar wind streams. 
We discuss here the three dimensional geometry of these two classes of interplanetary 
shocks and how these geometries can be studied with an out-of-the-ecliptic mission. 
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Out-of-the-EclipHc Distortion of Solar Blast Waves 

Lack of !.pherical s/mmetry in the solar wind distorts the surface of a shock 
wave as it propagctes away from its solar origin into interplanetary space. This 
phenomenon is observed near the solar equatorial plane at 1 AU, where the inhomo- 
geneities associated with solar wind streams produce typically a 30° deflection of 
the shock wave normal away from the radial direction ( Heinemann and S lscoe, ]974; 
Hirshberg et o l., 1974). The distortion results from differential advecrion of the 
shock front due to solar wind speed variations and from differential propagation 
speed of the shock in the solar wind due to density variations — the shock propagates 
more slowly in ' gh density regions. In spite of the large distortion of individual 
shocks, the averaged shock normal direction near the equatorial plane at 1 AU is 
radial from the sun ( Chao and Lepping , 1974; H einemann and Siscoe , 1974). 

A systematic variation In latitude oi fhe solar wind speed and density pro- 
duces a systematic distortion in latitude of the surfaces of solar produced shock waves. 
That is, the averaged shock normal direction (n ) in general will not be radial from 
the sun. The angle between ^ and the radial direction will depend on latitude in 
a manner which reflects the average latitudinal dependence of the solar wind speed 
and density. 

A lower iim t on the deviation of n^ from radial is shown in Figure 1. Here 
a shock wave that was spherical at 20 solar radii becomes distorted into a quasi- ellipse 
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at 1 AU by the action of a differential advection of 400 icm sec ct the equator 
increasing smoothly to 600 km sec ^ at the poles. This latitudinal gradient of solar 
wind speed, approximately 2 km sec ^dag is at the lower end of the range of 
gradients suggested in the literature (see reviews by Gosling , 1975, Hundhausen, 1975, 
Dobrowolfiy ond Moreno , 1975). The figure probably also under esti motes the dis- 
tortion for the assumed gradient since it neglects the possible latitudinal gradient in 
density-decreasing toward the poles (Hundhausen et al., 1971) — causing the shock to 
propagate faster in the solar wind at higher latitudes. 

The maximum deviation of n from radial (A6 ) occurs at mid-latitude a-d 

s max 

is about 25° at 1 AU. In order to observed this effect, enough solar generated shock 

waves must be measured while a spacecraft is in dte mid^latitude region to obtain a 

value for the polar angle (d^) of with a statistical error considercbly less than 

Ad . For the soke of having a numerical e.'.ampte, we take the requirement that 
max 

the expected standard deviation of 6^ be less riian or equal to 1/3 of the conservative 

value obtained above for A6 # i.e. S. D. (Q ) s8°. 

max s 

Near the equo^orial plane the distortions produced by solar wind streams 
cause approximately a 30° standard deviation in the or>gle between individual shocks 
and their average, radial direction. The standard deviotion of the average angle of 
N siiocks is 30 //N . For N = 14, this is ~ 8 . Thus In a 30 backgrounc 'oise 
level ir the angle of individual shocks due to stream structurevit takes 14 shocks to 
obi-^ir a value of 6 with an expected standard deviation o^ 8° from the true value. 

It is likely that the noise level at mid-latitudes •<! less than 30° since the effect of streams 
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in producing soior wind inhomogeneities as measured by enhanced radio sc inH I lotions 
is apparently confined to within 40*^ of the equator (Houminer^ 1973). Thus the 
requirement of 14 shocks is prob(4>ly more than necessary. 

To estimate the time needed to observe 14 solar generated shodcs, we use 
the 100 years of SSC stotistics compiled by Mayoud ( 1975 ), Xhe validity of this 
procedure is based on the study of Chop ond Lepping (1974) showing that at least 
87^ of SSC's can be associated with solar activity such as solar flores and type 2 
and type 4 radio bursts. There ore on average about 10 SSC's per year during solar 
minimum and about 35 SSC's per year during solar maximum. Thus about 17 months 
ore required in the first case and 5 months in the second of mid-latitude observations 
to observe 14 ^ock waves originating from a solar surface source. We conclude that 
an out-of-the-ecliptic mission scheduled for solar maximum would have high probd>ility 
of observing systematic shock wave distortion even if the worst-case example discussed 
above should apply. 

The importance of measuring the systematic distortion of the shock shape lies 
in its use in determining the systematic latitudinal dependence of solar wind parameters. 
This method is independent of all other methods. It does not involve unravelling 
separate space and time variations. In the case of on out-of-rhe-ediptic mission vio 
Jupiter, the space dependence involves both radial distance and latitude angle which 
ca»' give independent contributions to any variation: observed in in situ measurements. 
Knowledge of the three dimensionol shape of blast waves is important also for determining 
the flow of flare energy into interplanetary space. A non-spherical shock shape implies 
that , crgy is not distributed uniformly but converges in some places — relative to 
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a purely radial flow-- artd diveiges in others. In the case considered, energy would 
diverge away from the equator «o\»ing the shodc strength to decrease faster in the 
equotoriol plane thun would be expected on the basis of spho-ical symmetry. 

Three Dimensional Structure of Corotating Shodcs 

Near the equatorial plane the border between contiguous solar wind streams 
is a spiral (Scydahai , 1963; Dossier , 1967). If the trailir^ stream — in the sense 
determined by die direction of the solar rotation — is faster, a pair of shock waves 
will form at some distance from the sun (HuncBiousen , 1973 a,b). Such shock wave 
pairs have a parently been observed between 1 and 5 AU by the Pioneers 10 artd 11 
spacecraft (Hundhousen and Goslir^ , 1975; Smith ond Wolfe , 1975). In a steady 
state situation, the stream*;, their spiral border, and the shock pair ail corotate with 
the sun. In this section we estimate the heliocentric distance of shock formation in 
the equatorial plane as a function of the speed differential between the streams and 
give a qualitative description of the three dimensional shape of the shock surfaces. 
For the latter we assume that the stream border is perpendicular to the equatorial 
plane. This example illustrates the essential aspects of the geometry and possibly 
represents the typical case as indicated by radio scintillation observations (Houminer , 
1973) and the north-south alignment of coronal holes which might be the sources of 
fast streams (Krieger et al . , 1973; Noyes , 1975). 

Figure 2 shows the relevant geometry in the equatorial plane. The figure 
also illustrates one argument for expecting the existence of shock pairs which at the 
same time suggests a simple calculation for the opproximate heliocentric distance to 
their point of formation. The spiral labeled sfreom interface is the border between 
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the streams. Sample flow streamlines'** labeled fast and slow*-are shown in the corotating 
reference frame in which all geometrical features ore time independent. With the fast 
stream trailir^, the pitches of the spiral streamlines in both streams ore such Hiat they 
would intersect the spiral interface unless prevented from doing so by forces diat act to 
deflect the flows away. The build up of pressure at the interface resulting from the con-* 
vergence of the flow there produce such a force ( Siscoe , 1972). Relative to rite flows 
in the two streams, the interface looks, like a curving wall and rite flows are forced to 
follow the curve because of the increased pressure. Ihe compressive deflection of a 
sup^onic flow by a curving wall is knovra to produce a detached shock wave in the 
flow ( Londou and Lifschitz , 1959, p. 429). Streamlines intersecting the shock waves 
ore deflected p;»allel to the interface spiral, bringing to an end the compressive 
interaction between the streams. 

The oppr<»cimate location of the origin of the shock waves can be found by 
considering the characteristics of the flow eminating from a point. A, on the interface 
close enough to the sun that the streamlines are essentially radial but far enough from 
the sun to be in the supersonic region. The choract eristics are generated by following 
the progress of a sound wave starting at A and subsequently expanding and being con- 
vected with the flow, as illustrated in the figure. However, before the shocks ore 
formed, the flow converges cn the interface; thus fhe sound speed must be greater than 
the speed characterizing the convergence in order for the sound wove to expand. As 
the wove moves out from point A, the sound speed decreases because the solar wind 
cools as it expands, and the speed of convergence increases because of the relative 
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pitches of the spirals. A distance is reached when convergence exceeds the speed of 
sound and the wave begins to shrihk. 

From the point of view of an observer mowng with the flow along a streamline, 
and looking at the wall represented by the intmface, he sees the wall approach him— > 
that is, convergence in his frame of reference. If we think of the wall os a piston 
moving into the flow, at the point where the sound wave begins to shrink, tl.? piston 
is moving faster than the local speed of sound, and a shock wave will form upstream 
from the piston. Thus, the origins of the shock waves will be approximately at points B 
in the slow stream and C in the fast stream marking where the sound wave stops 
expanding owoy from the interface. 

To find the opproximote locations of these points and t'l. ir dependence on the 
speed oTfference between the streams, we consider an idealized case in which the slow 
stream has constant speed V^, the fast stream has constant speed Vp the speed at the 
interface is = (V^+Vf)/2, and the Mach number, M, is constant throughout. Using 
the procedure given in Heinemonn and Siscoe (1974), we find the equations of the 
characteristics to be 


n^=^(1+fcr/r^-r/r^) 
t?f = ^0 -0nr/rj-f/f^) 


0 ) 


where rj is the azimuthal angle in the corotating reference frame, r is the heliocentric 
distance, and = V r^ = V^/Mn . Without the in terms, these are the equations 
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of Porker's solar wind soirols ( Porker ^ 1963, p. 138). The 0n terms represent the 
movement of the sound wove away from the spiral. The equation of the interfoce is 

with ~ V n . To estimate where the ^ocks form we determine where the radial 
separation between the characteristics and the interface begins to decrease, i.e. set 
dAr/di 7 = 0. The result is shovm in Table 1 for different speed differentials, Vj-V^ 

_i 

with V = 400 km sec and M = 4. 
o 

As expected, the bigger the differential, the nearer the sun the shocks form. 

Ill the biggest cose considered 160 km sec \ they form near the orbit of Mors. Any 
differential bigger than approximately 50 km sec ^ produces shocks inside the orbit 
of Jupiter. Although this example is idealized, it gives a fair test of the argument 
based on chmocteristics to predict qualitatively the essential geometrical aspects 
of the formation of corotating shocks. The prediction that typical solar wind streams, 
which have differenticls bigger than 50 km sec should form shocks between the 
orbits of Earth and Jupiter is apparently confirmed by the Pioneer 10 and 11 observations. 
This justifies applying the argument to determine Hie out-of-the-ecllptic shape of 
cororcting shock waves. 

The application is straightforeward, and the result is immediate if we consider 
the situation at the poles. Here the border between thn streams is a radial line — the 
polar axis. The pitches of the streams lines are essentially zero and, hence, so is the 
speed of convergence of the streams. A sound wave starting here will expand forever. 
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although it wifi slow down because of the radiol decrease in sound speed. Ifius 
corotating ^odcs will not form over the poles. At intermediate latitudes, the pitches 
of the streamlines and the border are not zexo, but they are less than at the equator. 

A sound wave must travel further before the convergence speed overtakes the sound 
speed and diocks form. Thus, the distance to the formation of the shodcs increoses 
with latitude. 

The three dimensional geometry of the stream interface and the shodc pair 
is sketched in Figure 3. The interface is genmoted by on expandir^ meridian circle 
that rotates about the polar axis os it expands so that its intersection with the equatorial 
plane moves along the border spiral. The leading points of the shocks are in the 
equatcxial plane and the leading edges spiral outward os they move away from the 
equator, maintaining a proximity to the interface. A sketch on on expanded scale 
of o single corotating shock surface is shown in Figure 4. The motion of this surface 
is analogous to that of a tapered paper banner attached to a stick that is being twirled 
around the polo- axis. The length of the banner is determined by the lifetime of the 
solar wind stream. 

To study these structures an out-of-the-ecliptic mission is needed that covers 
the radial range between Earth and Jupiter. Knowledge of the three dimensional 
nature of th^e shock waves is essential for the interpretation of cosmic ray data and 
for applications to other astrophysicai situations. The distortion of solar shock wavK 
as described earlier and the highly structured geometry of corotating shocked woves 
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illustrate the complexity of the problems faced by galactic cosmic rays as they try to 
enter the inner solar system. Three dimensional probing of the interplanetary medium 
is required to obtain a complete picture of the extended stellar envelope of a 
representative from a major population of main sequence stars. The interaction of 
such stars with the interstellar medium in various galactic situations can they be 
treated with a fuller understanding of the stellar parameters. Comprehension of the 
three dimensional aspects of structures generated because a star rotates has application 
to contemporary ostrophysical problems such as the interaction of the Crds pulsar with 
the Crab nebula. 
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V--V 
f s 

dAr/dt} = 0 

16 (km/sec) 

12.5 (AU) 

32 

6.2 

48 

4.2 

64 

3.1 

80 

2.5 

96 

2.1 

112 

1.8 

128 

1.6 

144 

1.4 

160 

1.3 


Table 1. Approximate heolcentric distances in the equatorial plane to the 
formation of shock pairs dcj to interacting streams with various 
speed differentials. In this example the average speed is 400 km sec 
and the Mach number is 4 throughout. 
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Figure Captions 


Figure 1. A meridian plane cross section through a solar origin shock surface assumed 

to be circular at 20 solar radii and distorted into a quasi-ellipse at 200 solar 

radii by the action of equator to pole solar wind speed differential of 200 km 
-1 

sec . 

Figure 2. A sketch of geometrical features of the flow in the equatorial plane ond in 
the corotating reference frame. The stream interface separates a slow stream 
(leading) and a fast stream (trailing). Sample streamlines ond corotating shocks 
in both streams are shown. The shocks deflect the streamlines parallel to the 
stream interface — with no deflection they would intersect the interface. The 
circles are sequential snapshots of a sound wave starting at A and expanding 
while being convected with the flow (distortion due to the flow speed dif- 
ferential is neglected). Shocks form where the sound wave begins to shrink 
which happens when the speed of convergence of the streams toward the 
interface (due to the differences in the pitches of the various spirals) becomes 
bigger than the local sound speed. 

Figure 3. A sketch of the three dimensional geometry in the northern hemisphere of the 
stream interface and associated shock pair. The interfoce is generated by an 
expanding meridian circle that rotates obout the polar axis as it expands so 
that its intersection with the equatorial plane follows the interface spiral, like 
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Q needle in the groove of a record. The interface separates the shock po!r. 
Their leading point is in the equatorial plane and their leading edges spiral 
owciy from the sun as they move away from the equatorial plane, maintaining 
a proximity to the interface. 

Figure 4. A sketch on on expanded scale of a single shock surface. The location of 
the equatorial plane is indicated by lines radiating from the sun. The motion 
of the surfoce is similar to that of a tappered paper banner attached at its 
point to a stick that is being twirled around the polar axis. The length of 
the banner is determined by the lifetime of the streams. 
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